Finite size effects on kaonic pasta structures 
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Non-uniform structures of mixed phases at the first-order phase transition to charged kaon con- 
densation are studied using a density functional theory within the relativistic mean field model. 
Including electric field effects and applying the Gibbs conditions in a proper way, we numerically 
determine density profiles of nucleons, electrons and condensed kaons. Importance of charge screen- 
ing effects is elucidated and thereby we show that the Maxwell construction is effectively justified. 
Surface effect is also studied to figure out its effect on the density profiles. 



I. INTRODUCTION 

Concerning phase transitions in nuclear matter, var- 
ious scenarios have been discussed, like e.g. liquid-gas 
transition pj , pion condensational, hadron-quark decon- 
finement transition 0, El HI El 01 HQ* transitions between 
different color superconductin g p hases [ToL ITU IT^ . fl3| . 
quark ferroma gne tic transition |14| , charged p meson con- 
densation [TEl lltH. etc., and kaon condensation. In most 
cases these are first-order phase transitions (FOPTs). 

References E3, EH |2lL 13 H3, H3| 

and many others studied possibilities of s-wave kaon 
condensations, whereas Refs. [2I considered also 
those of p-wave condensations. It has been argued (see 
ElHEIillilllillilliiM and refer- 

ences therein) that the FOPT to a K~ condensate state 
might occur in neutron stars already at densities only 
several times larger than the nuclear saturation density 
Po- 

The possibility of the FOPT to the kaon condensate 
state may lead to interesting consequences for physics of 
neutron stars, such as a delayed collapse of a supernova 
to a low-mass black hole, a fast cooling mechanism of 
neutron stars due to the opening up of the nucleon direct 
Urea process under a background of the charged kaon 
condensate, etc [13, IS3, ISH E3, E3, EE IS^ - 

Reference argued that in systems with different 
charged species the structured mixed phases might ap- 
pear due to the FOPT. For example in case of the nuclear 
matter below the saturation density, where the liquid-gas 
FOPT is relevant, the "nuclearpasta" structures have 
been studied by many authors [1 El El El El El El 
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The physical reason for the possibility of the mixed 
(pasta) phases with different geometrical structures is 
that the charge neutrality may hold only globally rather 
than locally as in the Maxwell construction. Mechani- 
cally, a balance between the Coulomb force and the sur- 
face tension is responsible for the occurrence of the spa- 
tially non-uniform structures. 

The appearance of the pasta structures could have im- 
portant consequences for the various neutron star phe- 
nomena. If appeared at an early stage of the neutron 
star evolution the pasta could cause a drastic change of 
the neutrino opacity. It may influence the subsequent 
neutron star cooling, it affects the matter resistance to 
the stress and consequently the glitch phenomena, etc. 

Long time there existed a naive view that not all the 
Gibbs conditions can be satisfied in a description by the 
Maxwell construction H l23l l24l l25j . because the local 
charge neutrality is implicitly assumed in it. As the re- 
sult of previous works, cf. Hi mil, it was suggested 
that a broad region of kaonic pasta as a mixed phase may 
occur in neutron stars with various geometrical struc- 
tures of the kaon condensed phase (dense phase) embed- 
ded in usual nuclear matterphase (dilute phase). How- 
ever, in recent papers [100 H |H OH] we have demon- 
strated that both treatments using the Gibbs conditions 
and the Maxwell construction are actually in agreement 
with each other, if one properly includes electric field ef- 
fects. The bare Coulomb field is screened by rearranged 
charged particles. Taking first the hadron-quark phase 
transition as an explicit example [1H0, we showed that 
with screening effects included the most realistic param- 
eter choice (density, chemical potentials, surface tension, 
etc.) does not allow the mixed phase due to the mechan- 
ical instability of structures. Thus if exists, the density 
region of the structured mixed phase should be largely 
limited. 

In subsequent papers H El Ell > we have explored the 
charge screening effect on the nuclear pasta phases at 
sub-nuclear densities. We compared the full calculation 
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with those either not including the screening at all or in- 
cluding it only non-systematically, as has been performed 
in previous works. We found that the screening effect, be- 
ing systematically treated, significantly affects the den- 
sity range of the nuclear pasta with various geometries, 
while it results in only minor corrections to the EOS. 
We have also emphasized that the proton spatial rear- 
rangement effect resulting in a screening is much more 
efficient, compared with the electron screening, because 
the proton Debye screening length is essentially smaller 
than the electron one. However, since the charge density 
as well as the baryon number density are rather low in 
the nuclear pasta case, the screening effect is still not so 
prominent compared to that we find for denser systems. 

In this paper we continue our research 
of the special role of the charge screening effects in the 
problem of the construction of the mixed phases, tak- 
ing s-wave kaon condensation as a further example, cf. 
brief notes [Hfll l5lj . We study whether the kaonic pasta 
appears in dense neutron star matter, and examine the 
validity of the Maxwell construction for the practical use. 

One may believe that models taking into account the 
chiral symmetry are more realistic for the treatment of 
the kaon condensation. Actually man y w orks have been 
done using this approach [IS l2fl 123 • However, 

we will exploit here a non-chiral model since the bulk 
calculations performed within chiral models do not allow 
the mixed phase pH I52L Previous investigations of 
the kaon mixed phasejUHl EI El E! El were per- 
formed in a framework of the RMF models. Therefore 
we will also exploit a RMF model for a more straight 
comparison with the previously obtained results. 

We also will figure out the surface effect. There are few 
works about this matter so far j2f|, EE 53] , and the finite 
size effects of the Coulomb interaction and the surface 
tension are still not elucidated. In Ref. [2{| authors dis- 
cussed the surface effect and derived the surface tension 
between the normal nuclear matter phase and the kaon 
condensed phase by considering two semi-infinite matters 
omitting the Coulomb interaction. They suggested that 
the kinetic terms consisting of space gradients of the a, 
u>, p mean-fields affect only slightly the density profiles 
and the resulting surface tension. Therefore they omitted 
these kinetic terms while they kept the kinetic term of the 
kaon field. This scale of the kaon field then determines 
the change of all meson fields and the nucleon fields, as 
well as the resulting surface tension, see their Fig. I. 
Then, they find the surface tension which is now deter- 
mined by this single length-scale. After that, for different 
geometries minimizing the surface plus the Coulomb en- 
ergy per volume in the structure size one may find the 
optimum structure size and the geometry. Note that the 
Coulomb energy density in their approach is given by the 
bare Coulomb potential decoupling with other fields. In 
reality this approximation does not hold. The scale of the 
change of the kaon field proves to be of the same order of 
magnitude as the scale of the change of the electric field, 
i.e. the Debye screening length. A variation of the electric 



field significantly affects the kaon distribution and vise 
versa. Therefore, we must discuss the surface effect by 
the kaon field and other meson mean- fields, separately: 
only short-scale fields can decouple with long-scale fields 
and thus the contribution of the former fields may be re- 
placed by the phenomenological surface-tension parame- 
ter, whereas long-scale fields should be treated explicitly. 

In Ref. [2(| the authors omitted kinetic terms of the 
mean a, u>, p fields (neglecting also the corresponding 
surface tension effect) and discussed the charge screen- 
ing effect. By numerical analysis they demonstrated how 
much the screening effect may affect the density profiles 
for the kaon droplet phase. 

We consistently incorporate the Coulomb field as the 
gauge field coupled with other fields of charged particles 
and we keep kinetic terms of all mean fields. We demon- 
strate that the cr, u, p fields are the short-scale fields. 
Thus we treat consistently the long-range effects, such as 
screening, and the short-range effects and we show their 
interplay for the kaonic pasta. This we believe gives a 
new insight into the finite size effects in the mixed phase. 
Also by artificial variation of the scale of the short-range 
fields we demonstrate effects that can be reduced to the 
surface tension. 

The paper is organized as follows. In Sect, [n] we in- 
troduce the thermodynamic potential and equations of 
motion for fields. In Sect. ITTT1 we numerically solve the 
Poisson equation and the mean field equations to find 
the electric potential profile and the density profiles. We 
also derive the EOS and the phase diagram for the kaonic 
pasta. In Sect. IIVI we discuss the screening effect and 
the surface effect on the kaonic pasta. Our conclusions 
are drown in Sect. E] Details of the perturbative treat- 
ment of the Coulomb effect, which we refer to as the "no 
Coulomb" approximation, are deferred to Appendix lAl 



II. THERMODYNAMIC POTENTIAL AND 
NUMERICAL PROCEDURE 

A. Thermodynamic potential 

Let us present our framework to study kaonic pasta 
structures. Following the idea of the density functional 
theory within the RMF model, we can formulate equa- 
tions of motion to study non-uniform nuclear matter nu- 
merically, cf . [54| . The RMF model with fields of mesons 
and baryons introduced in a Lorentz-invariant way is rel- 
atively simple for numerical calculations, but on the other 
hand it is sufficiently realistic to reproduce bulk proper- 
ties of finite nuclei as well as the saturation properties of 
nuclear matter [4^ . In our framework, the Coulomb in- 
teraction is properly included in equations of motion for 
nucleons and electrons and for meson mean fields, and 
we solve the Poisson equation for the Coulomb potential 
Vboui self-consistently with those equations. Thus the 
baryon and electron density profiles, as well as the me- 
son mean fields, are determined in a way fully consistent 
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with the Coulomb interaction. 

We start with the thermodynamic potential for the sys- 
tem of neutrons, protons, electrons and mesons including 
kaons 



(1) 



The first term 



n N = J2 



d 3 7 



fk 
47T 3 



N 



k 2 



PaV a 



(2) 

is the contribution of nucleons with the local Fermi mo- 
menta kp a (r);a = n,p, m* N (r) = m^r — g<jN&{ Y ) is the 
effective nucleon mass and hin is the nucleon mass in the 
vacuum. Nucleons couple with a, lo and p mesons and 
thereby, 

v n (y) = p n - ,g w jvw (r) + g pN R a (r), (3) 
v p (t) = fi p + Vcoui(r) - guN^o(r) - g pN R (r), 

where p n and p p are neutron and proton chemical po- 
tentials and ga-N, g^N and g p x are coupling constants 
between corresponding fields. 

The second term in incorporates the scalar (a) and 
vector (u>q,Rq) mean fields, 
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d 3 r 



(Va) 2 



(Vwo) 2 + mlu; 2 



- + U(a) 
(Vi? ) 2 + m 2 R 2 
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where m CT , m u and m p are the field masses, and U(a) = 
■^bm]s[(g a NO') 3 + jc(g<rNO') 4 is the nonlinear potential for 
the scalar field. 

The third term in Q contains the contribution of the 
Coulomb field (described by the potential Vc ou i( r )) and 
the contribution of relativistic electrons, 



d A r 
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where p e is the electron chemical potential. 

The last term in is the thermodynamic potential of 
the mean K~ meson field, 



9. 



K 



d A r 



2 

gu)K^>o 



' K 



g P KRo) 



+ (pK — Vcoul 



/i(V0) 2 



(G) 



where m* K (r) = mx — 5o-/fc(r) is the effective K~ mass, 
ttik is the kaon mass in vacuum, g a K, g^K, gpK are 
the coupling constants, px is the K~ chemical poten- 
tial, K = f K 0/V2 is the kaon field and f K ~ 93 MeV is 
the kaon decay constant |69|. The kaon charge density 
Pk is expressed in terms of the kaon field as 



PK 



(p K - Vboui + guK^o + g P KRo) Ik® ■ ( 7 ) 



Temperature T is kept zero in the present study. 

For nucleons and electrons we used the local-density 
approximation, i.e., nucleons and electrons are described 
by their local densities. This approximation has its sense 
only if the typical length of the change of the nucleon 
density is larger than the inter-nucleon distance. Deriva- 
tive terms of the particle densities can be incorporated 
in the quasi-classical manner by the derivative expansion 
within the density functional theory Their contri- 

bution to the energy can be reduced to a surface tension 
term. Here we simply discard those derivative terms, as a 
first-step approximation. Thus we discard the contribu- 
tion of the nucleon fields to the surface tension assuming 
that it is smaller than the corresponding contribution of 
the meson fields that we retain. In the case when we sup- 
press derivative terms of the nucleon densities they follow 
changes of the meson a, lo, p mean fields and the kaon 
and the Coulomb fields that have derivative terms. Note 
that we have fitted our model to properly describe finite 
nuclei (see below) without including nucleon derivative 
terms. If we introduced them it would need to re-adjust 
the model parameters. We should bear in mind that for 
small structure sizes quantum effects become prominent. 
For simplicity we disregard these effects. Thus within 
this scheme we may properly describe only rather large- 
size structures. 



B. Equations of motion and numerical procedure 

Equations of motion for the mean fields including kaon 
and for the Coulomb potential are obtained from the vari- 



ational principle: 



sn 



= 



: <r,UQ,Ro,0,Vcoui) and 



Sj^rpj =0 (a = n,p, e). These equations read 
V 2 (i(r) 



+ 2 g<7K m K f 2 K 2 (r), (8) 
V 2 w (r) = mluj (r) - g^Nippir) + p n (r j) 

- f K guK0 2 {r){pK - Vcoui(r) 
+9wkUq{*)+9 p kRo{*)), (9) 

V 2 i? (r) = m 2 i? (r) - g pN (p p (r) - p n (r)) 

- fKg P Kd 2 (r)(p K - Vboul(r) 

+ ffwJfWo(r) + g pK Ro{r)), (10) 
V 2 0(r) = [m* K (r) 2 - {p K - V Co ui{r) 

+ 9uKUa(r) +g pK Ro{r)) 2 ]0{r), (11) 

W Coul (r) = 47re 2 p ch (r), (12) 
the charge density p c h(r) = p p (r) + p e (r) + p K {r), 

Pn = Pb = \j k Fn (r) 2 + m* N (r) 2 

+ .g^A^o(r) - </ P Ari?o(r), (13) 

- Pe = \J k Fp (r) 2 + m* N (r) 2 
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TABLE I: Parameter set used in RMF in our calculation. 



ffo-iv ff^iv ffpjv b c rna rn^j m p fn) rn K g^K g P K Uk(po) 



6.3935 8.7207 4.2696 0.008659 0.002421 400 MeV 783 MeV 769 MeV 93 MeV 494 MeV g^ N /3 g pN -120 - -130 MeV 
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FIG. 1: The mass number A dependence of the binding energy of finite nuclei (left) and the density profile of 208 Pb nucleus 
(right). Different curves correspond to different meson masses used in the RMF calculation. Thick horizontal gray line in the 
left panel shows the case of the nuclear matter for the comparison. 



+ g^NOJ^v) + g pN R (r) - Vcoul(r). 

(14) 

The last two equations are the standard relations between 
the local nucleon densities and chemical potentials. We 
have assumed that the system is in the chemical equilib- 
rium in respect to the weak, electromagnetic and strong 
interactions and we introduced the baryon chemical po- 
tential /is = fi>n an( i the charge chemical potential, i.e. 
the electron chemical potential, /i e , according to the cor- 
responding conserved charges. Under the same assump- 
tion fJ, K = (J, e - 

To solve the above coupled equations numerically, the 
whole space is divided into equivalent Wigner-Seitz cells. 
The geometrical shape of the cell changes as follows: 
sphere in three-dimensional (3D) calculation, cylinder 
in 2D and slab in ID, respectively. Each cell is glob- 
ally charge-neutral and all physical quantities in the cell 
are smoothly connected to those of the next cell with 
zero gradients at the boundary. Every point inside the 
cell is represented by the grid points (number of grid 
points iVgrjd ~ 100) and differential equations for fields 
are solved by the relaxation method for a given baryon- 
numbcr density under constraints of the global charge 
neutrality. Details of the numerical procedure are ex- 
plained in Ref. |49| . 



C. Parameter set and finite nuclei 

For the study of a non-uniform nuclear matter, the 
ability to reproduce the bulk properties of finite nuclei 
should be essential. Parameters of the RMF model are 



chosen to reproduce saturation properties of nuclear mat- 
ter: the minimum energy per nucleon —16.3 MeV at 
p = po = 0.153 fm~ 3 , the incompressibility K(po) = 240 
MeV, the effective nucleon mass m* N (po) = 0.78mAr; 
tojv = 938 MeV, and the isospin-asymmetry coefficient 
Osym = 32.5 MeV. Coupling constants and meson masses 
used in our calculation are listed in Table I. 

The parameter g a K enters the value of the K~ optical 
potential Uk defined by Uk — 9oK<3 + 9ujK^o- There 
have been many works trying to fix Uk at the saturation 
density from the data on the kaonic atoms [5^, |5(j and 
from calculations but there is still a 

controversy in its depth. We take here a somewhat deep 
potential, as shown in Table I, to compare our results 
with the previous ones [U 0, H^, |2(J . To understand 
a dependence of the results on the value Uk we further 
allow for a variation of it. 

We have checked that with the meson masses listed 
in Table I (m a = 400 MeV, etc.) and by including the 
Coulomb interaction, the binding energies of finite nuclei 
and the proton fraction, as well as the nucleon density 
profiles are well reproduced, except for very light nuclei 

We examine in this section how the surface tension 
could be correctly incorporated in our RMF calculation. 
We explicitly treat the kaon gradient term in (|llf> since, 
as we show, the kaon and Coulomb fields are essentially 
coupled and due to that the scales of changes of these 
fields prove to be of the same order of magnitude. These 
are long scales. Reference using an analytical model 
describing the charge distribution in a pion condensate 
droplet has demonstrated that actually there exists only 
one long scale in the problem. Namely due to that one 
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FIG. 2: (Color online) Density profiles of kaonic structures. 
Here the density does not mean charge-density but number- 
density of particles. Uk, the kaon optical potential at the 
nuclear saturation density, is set to be —130 MeV. Displayed 
using the right axis is the Coulomb potential Vcoul (written 
as "V"). 



cannot reduce the kaon contribution to a purely surface 
term in case the Coulomb interaction is introduced in 
equations of motion. 

On the other hand the a, uj, p prove to be short-range 
fields. In order to test effects of gradient terms of the 
short-range meson fields in Q and (|S|)- (|10fl we multiply 
meson masses m a , and m p by a factor cm- We take 
cm = 1 (in realistic case), 2.5 and 5.0. For example, 
with the factor cm = 2.5 we get the choice m (T =1000 
MeV, m w =1958 MeV and 77^=1923 MeV. However to 
obtain appropriate matter properties, we simultaneously 



fix the ratio gi N /i 



-a , uj and p) for all values of 



because thermodynamic characteristics depend only on 
9%nI w ^> rather than on g^ and separately. 

First to show the effect of thus simulated surface ten- 
sion we discard the Coulomb interaction (except that we 
use the global charge neutrality condition). Left panel 
of Fig. n demonstrates the binding energy of finite nuclei 
calculated with different meson masses. We have checked 
that for all sets of a, uj, p meson masses (for different cm) 
the curves approach the value of the binding energy of 
the nuclear matter for A — > oo. Also by the use of heavy 
meson masses, the binding energy of finite nuclei (for fi- 
nite A) approaches to that of nuclear matter indicated by 



a thick gray line. This shows that the surface tension is 
reduced with the increase of the meson masses, cf. [63| . 
Notice that this statement is correct only if we fixed the 
ratio g\ N lm\. 

Then we allow for the Coulomb interaction. Right 
panel of Fig. ^ shows the baryon density profile of 208 Pb 
nucleus obtained with different cm factors. With a larger 
cm, i.e. with a weaker surface tension, the baryon den- 
sity at the center decreases and that near the surface 
increases due to the Coulomb repulsion of protons. 

When we further discuss the kaonic pasta structures in 
the later section, we compare results obtained with stan- 
dard meson masses and those with heavy meson masses 
to figure out the effect of the surface tension due to the 
meson mean-fields. 



III. RESULTS 

The density profile is determined by solving equations 
of motion (|S|)- (|14[) . For each baryon density, the cell ra- 
dius and the geometrical dimension are chosen to mini- 
mize the energy. If Glendenning's claim were correct, the 
structured mixed phase would develop in a broad density 
range from well below to well above the critical densities 
determined by the Maxwell construction. In this density 
interval the matter should exhibit the structure change, 
similar to the nuclear "pasta" phases jl^: the kaonic 
droplet, rod, slab, tube, bubble. We observe such struc- 
tures in our calculation as well. 

Figure [21 displays typical density profiles and the 
Coulomb potential. The horizontal axis is the distance 
from the cell center and the hatch shows the cell bound- 
ary. The symbols "3D" (three dimensional) etc. indi- 
cate the dimensionality of the geometry. From the top 
of the figure, the matter structures correspond to kaonic 
droplet, rod, slab, tube, and bubble. The neutron distri- 
bution proves to be rather flat. The Coulomb potential 
is given by imposing the gauge condition. Here we fix the 
gauge by the condition, Vcoui(-Rccii) = for all the cases. 
Note that the maximum value of the Coulomb potential 
becomes rather large. 

In the upper panel of Fig. [3] we depict the energy per 
nucleon of the matter. The dotted line indicates the 
case of single phase (if one assumes absence of the mixed 
phase). In this case the uniform matter consists of nor- 
mal nuclear matter for undercritical densities and the 
kaonic matter for overcritical densities. The cross on the 
dotted line (pb — 0.46 fm -3 ) shows the critical density, 
i.e. the point where kaons begin to condensate in the 
case of single phase. Pieces of solid curves, on the other 
hand, indicate energetically favored structures. Droplets 
begin to appear for ps > 0.41 fm -3 smoothly decreas- 
ing the energy of the system. The mixed phase disap- 
pears for pb > 0.74 fm -3 . The occurrence of the kaonic 
pasta structures, i.e. kaonic droplet, rod, slab, tube, and 
bubble, results in a softening of the matter (the energy 
decreases with their appearance). 
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0.4 0.5 0.6 0.7 0.8 
Pb [fm~ 3 ] 

FIG. 3: (Color online) Top: Binding energy per nucleon of 
the nuclear matter in the beta equilibrium. The dotted line 
below the cross shows the uniform normal nuclear matter and 
above the cross, the uniform kaonic matter. Bottom: Struc- 
ture size Rd (thick curves below) and cell size 7? ce ii (thick 
curves above). Compared are the Debye screening lengths of 
electron, proton and kaon. The Debye screening lengths are 
calculated using the explicit dependence of p a on p a ■ 



hand, are always shorter than i? C eii and the structure size 
Rd- Now we are able to estimate the value of the result- 
ing total Debye screening length A^j *"*. For example in 
the case of ps = 0.55 fm~ 3 we obtain A^j *-* ~ 3 fm. 

On the other hand the typical length of the Coulomb 
potential is ~ 5 fm for ps = 0.55 fm -3 (see Fig. 01. This 
value is longer than the above estimated Debye screening 
length due to the implicit dependence of p a on p\, a - 

From this comparison we conclude that a consistent 
inclusion of the Coulomb screening effect is not a triv- 
ial problem and it has really important consequences for 
the correct description of density profiles of protons and 
kaons (see Fig. below). 

So far we have presented results for the K~ optical po- 
tential in nucleus equal to Uk = —130 MeV. With a lower 
value Uk = —120 MeV, the density range of the kaonic 
pasta is narrowed: from ps — 0.49 to 0.71 fm -3 instead 
of from pb = 0.41 to 0.74 fm -3 in the case of Uk = —130 
MeV. The effect of the kaon pasta structures on the EOS 
(the energy difference between single phase calculation 
and full calculation) becomes still smaller. The kaon and 
proton densities in a droplet become 20 - 30 % lower than 
those for the case Uk = — 130 MeV, and the neutron den- 
sity increases. Other features remain qualitatively the 
same. 



IV. CHARGE SCREENING AND SURFACE 
EFFECTS 



In the lower panel of Fig.|3|we plot the structure sizes 
(for droplet, rod, slab, tube and bubble) Rd and the cor- 
responding cell sizes i? C cii- The structure size Rd is cal- 
culated by the use of the averaged density and its spacial 
fluctuation as in Ref. |4^. We find that at the onset 
density the radius of the cell is infinitely large in case of 
the full calculation. The corresponding steep increase of 
■Rceil with decreasing density is clearly seen in the figure. 

To roughly estimate the Debye screening length Ap 
(see Appendix [E] for detail) let us for simplicity use 
only explicit dependence of (p a ) on p a in (|B3|) . i.e., 



let us replace < —^£®l by 



dp a 



dp a 



hence 



d (p p ) ( k F P ) \l (k Fp ) 2 +m N 2 



dp e 



dp p 



for electron and proton, 



and 



d(pK) 
dp, 



K 



{—f 2 8 2 ) for kaon. Here, the bracket 



indicates an averaging over an appropriate region (cf. Ap- 
pendix [Bj. The dashed lines and the dotted line in the 
lower panel of Fig. show thus treated partial contribu- 
tions to the Debye screening lengths of the electron, pro- 



ton and kaon, A^, A^, and A^ ; , respectively. We see 

that in most cases A^ is less than the cell size i? C cii but 
it is larger than the structure size Rd- The proton Debye 
length A^ and the kaon Debye length A^ , on the other 



(if) 



A. Charge screening effect 

To demonstrate the charge screening effects we com- 
pare our results with those given by a "perturbative" 
treatment of the Coulomb interaction often used in the 
literature, "no Coulomb" calculation (see AppcndixlAlfor 
details). The electric potential is discarded in equations 
of motion l(7|l- lfTTl) . (JT3J, (|T4*jl which determine the den- 
sity profiles. The Coulomb energy is then added to the 
total energy by using the charge density profile thus de- 
termined to find the optimal value with respect to the cell 
size -Rcdi- Comparing Figs.|21and we see that the den- 
sity range of the mixed phase is more narrow in the case 
of the "no Coulomb" calculation than in the full calcula- 
tion, while the energy gain is almost the same. A remark- 
able difference is seen in the cell radii, especially near the 
onset density of kaonic pastas, for ps < 0.5 fm -3 . The 
cell size given by the full calculation is always larger than 
that given by the "no Coulomb" calculation. 

To elucidate the screening effect, we depict the i? C eii 
dependence of the energy per nucleon in Fig. [S] in both 
cases. In the full calculation a large cell with a small 
volume fraction / = (Rd/Rcc\i) 3 appears near the onset 
density, which situation is close to the uniform single 
phase. On the other hand, a smaller cell appears near 
the onset density in the "no Coulomb" calculation. The 
energy gain is higher in the full calculation. The cell 



7 



160 
> 140 
^ 120 

I ioo 

80 



S- 20 



10 




U K =-130 MeV 
no Coulomb 





droplet 
rod 
slab 
tube 
bubble 
uniform 
single phase 

I i i i i 



0.4 0.5 0.6 0.7 0.8 
Pb [fm~ 3 ] 



96 
94 
92 



I no 



i 105 



130 



125 



0.45 fm" 3 droplet 



— * — no Coulomb 
— * — full calc 




10 

flcell [ fm ] 



20 



FIG. 4: (Color online) Top: the binding energy per nucleon 
of the nuclear matter in beta equilibrium. The electric po- 
tential is discarded determining the density profile and added 
evaluating the energy. Bottom: the structure size Rd (thick 
curves below) and the cell size R ca n (thick curves above). 



radii obtained in both calculations deviate essentially for 
droplets, less for rods, still less for slabs, etc. 

The differences of both calculations are also seen from 
Fig. EJ where we present the pressure as a function of 
the baryon-number density. In the case of the full cal- 
culation the pressure changes continuously. In the case 
of the "no Coulomb" calculation there arises a jump in 
the pressure at the onset density at which the droplets 
begin to appear. The pressure jump clearly shows that a 
large structure change occurs at the onset density. The 
discontinuity arises as an artifact of the "no Coulomb" 
calculation. The Coulomb energy is overestimated in this 
case for large values of i? cc ii, compared to that for the 
full calculation, where occurs the screening effect. In- 
deed, the Coulomb energy per droplet grows proportion- 
ally to R^ cii in the case of the "no Coulomb" calculation 
since the bare Coulomb interaction has an infinite range. 
On the other hand in the full calculation there appears 
another length scale, the Debye screening length, and 
thereby the screened Coulomb interaction is no longer 
scale-invariant. Thus when the structure radius is sig- 
nificantly larger than the minimal screening length the 
Coulomb energy contribution becomes ineffective. 

In Fig. El we also depict the pressure when the Gibbs 
conditions are applied for two semi-infinite matters disre- 
garding the Coulomb interaction (indicated by "Gibbs" ) 
and that given by the Maxwell construction (indicated 
by "Maxwell"). We see that the pieces of solid curves 
lie between "Gibbs" and "Maxwell". The full calcula- 



FIG. 5: (Color online) The cell size f? C eii dependence of the 
energy per nucleon. Crosses on the curves indicate the mini- 
mum points. 



tion case is more similar to the one given by the Maxwell 
construction. 

To further demonstrate the charge screening effect on 
the kaonic pasta, We compare the density profiles ob- 
tained by the full and "no Coulomb" calculations in 
Fig. In case of the full calculation the difference be- 
tween the negative charge density (of kaons and elec- 
trons) and the positive charge density of protons is 
smaller, demonstrating that the profiles are more close to 
those given by the local charge neutrality. These results 
also suggest that the Maxwell construction is effectively 
justified in the full calculation case owing to the charge 
screening effects. Also we see that the absolute value 
of the kaon charge density is substantially larger in the 
"no Coulomb" case. This gives us an additional argu- 
ment for the essential coupling of the kaon and Poisson 
equations. Modification of the Coulomb term in the kaon 
equation significantly affects the charged kaon distribu- 
tion whereas changes in the proton, neutron and electron 
distributions are more smooth. 



B. Surface effects 

Next we investigate the effect of the surface tension 
caused by the meson mean-fields on the kaonic pasta 
structures. The previous calculation by Norsen and 
Reddy [2(1 153| was done discarding the gradient terms 
of a, lu and p mesons. The calculation with infinitely 
heavy meson masses would correspond to it. 

Figure shows the same quantities as Fig. O EOS and 
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length scales. The Coulomb term is mainly balanced by 
the kaon term, and the structure size is actually deter- 
mined by the minimum of the Coulomb plus kaon kinetic 
energy. The pasta structure is realized by the balance 
between the surface tension and the Coulomb repulsion. 
In this sense, tr, lo and p mesons have less contribution 
to the surface tension compared to kaon in the case of 
kaonic pasta structure, although the clear distinction of 
the Coulomb repulsion and the surface tension is diffi- 
cult due to the screening effects and the large surface 
diffuseness. 
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FIG. 8: (Color online) Top: Binding energy per nucleon of the 
nuclear matter in beta equilibrium. Bottom: Structure size 
Rd (thick curves below) and the cell size i? C cii (thick curves 
above). 



FIG. 7: (Color online) Comparison of density profiles of 
kaonic matter in full calculation and in "no Coulomb" case. 
The cell size, R cc n = 17 fm, is not optimized since the opti- 
mum values would be different for different treatments of the 
Coulomb interaction. 



the structure size, but now for the case of an artificially 
suppressed surface tension (for cm — 5.0). Comparing 
Figs. |H1 and |3 we see that there is almost no difference in 
the EOS. However, the density range of pasta structure is 
slightly broader for the case of the weaker surface tension. 
The reason for a similarity of Figs. |H1 and El is as follows: 
Even without the contribution of gradient terms of a, 
lj and p mesons (that corresponds to infinitely massive 
mesons, however for fixed g^ N / 'm|) , there is still a con- 
tribution from the kaon gradient term in Eq. (jl^). In our 
case the kaon and the Coulomb fields have similar typical 



V. SUMMARY AND CONCLUDING REMARKS 

We have discussed the effects of the charge screening 
and the surface tension on the kaonic pasta. Since the 
kaon mixed-phase appears at high-densities, we see that 
changes are more remarkable than for the "nuclear pasta" 
at subnuclear densities .4^. We found that the density 
range of the structured mixed phase is largely limited 
by the charge screening and thereby the phase diagram 
becomes similar to that given by the Maxwell construc- 
tion. Although the importance of such a treatment has 
been demonstrated for the hadron-quark matter transi- 
tion 0, , one of our new findings here is that we can 
figure out the role of the charge screening effect without 
introducing any "artificial" input of the surface-tension 
parameter. For this aim in our study we have used the 
relativistic mean field approach. 
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On the other hand we remind that the bulk calculation, 
where finite size effects (the Coulomb interaction and the 
surface tension) are not included, does not produce the 
mixed phase in the chiral models 0, I27I IH^ . Thus it 
should be interesting to perform the same self-consistent 
study with properly included finite size effects, as we have 
done above, but within a chiral model. 

In application to the newly formed neutron stars in su- 
pernova explosions, finite temperature and neutrino trap- 
ping effects become important, as well as the dynamics 
of the first order phase transition with formation of the 
structures. It would be interesting to extend our frame- 
work to include these effects. 

We focused our study on the s-wave kaon condensa- 
tion. However it would be interesting to extend it to the 
possibility of the p-wave kaon condensation. 

It might be also interesting to a ppl y our framework to 
the problem of kaonic nuclei |64l l65l lolj l67| , similar to 
the description of finite nuclei. Within our framework we 
can construct a consistent description of kaonic nuclei by 
taking into account the Coulomb interaction as well as 
the KN interaction 681. 
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APPENDIX A: "NO COULOMB" CALCULATION 

We briefly describe here the "no Coulomb" approxi- 
mation, which has been frequently used in the bulk cal- 
culations. To begin with we present the thermodynamic 
potential as 



^(Vcoul) = ^matter (Vboul) — 



1 



87re- 



■(vy Co ui) 5 



(Al) 

where we isolated the electric field contribution and 
^matter (Vboui) summarizes other baryon, electron, kaon 
and mean-field contributions. Expanding f2(Vboui) 
around a reference value Vcoui — Vo — const, we have 



o(v< 



CoulJ 



^matter (Vo) 
^matter (Vcoul) 



SV Coul (r) 



(Vboui(r) - Vo) 



87re 2 



+ O[(V C oui(r)-V- ) 2 ].(A2) 



The non-zero value of Vq shifts the value of the charge 
chemical potential in such a way that only the combina- 
tion u — Vboui + Vq is meaningful due to the gauge in- 
variance. To be specific in the paper body we performed 
calculations in the "no Coulomb" case with Vo = 0. Since 
^matter (Vboui) includes Vboul in a gauge invariant way, it 
can be further rewritten as 

fi(Vboul) = Onatter(Vb)- / d 3 rp ch (V ) (Vboul (f) - V ) 



d 3 r (VVbo U i) +0 ((vb oul (r)-Vb) 2 ). (A3) 
8ne z 



Using the total charge neutrality condition, it eventually 
reads 

n(Vbanl) = ^ NC (Vboui) + O ((VboulM - V ) 2 ) , (A4) 
with 



tt NC (Vb ou i) = n matt or(Vb) - / rf 3 r /9 ch (Vb)Vc ou i(r) 



d 3 r (VVboul) 



8ne 2 



(A5) 



The Poisson equation for Vboul is given by 

V 2 Vboui = 47r e y h (V ), (A6) 

while the equations of motion for other fields and density 
profiles, / Q (r) = <^(r), p a (r), read 



ffl ma ttcr(Vb) dp ch (V Q ) 
Sf a (r) df a (r) 



U Cou i (r) = 0. (A7) 



In (|A7|) one can neglect the second term of the electro- 
magnetic origin compared to the first term related to the 
strong interaction. Then one has 



^matter (Vb) 



0. 



Sf a (r) 

Using Eq. (|A6|) . we can recast Eq. (| A5|) into the form 

1 



(A8) 



(Vboul) = OnattorCVb) - g / r p ca {V )V Cou y(r) 



= ftmatter(Vb) + / <Tr 



,3.(VVboulT 

87re 2 



(A9) 



Equations 1A6|I . (|A8|) and l|A9|l are used in the "no 
Coulomb" calculation. The second term in Eq. l|A9Jl (let's 
call it the Coulomb energy) gives a positive (repulsive) 
contribution. 

The corresponding total energy E NC is then given by 
the Legendre transformation from il NC , 



Vo 



£NC =n NC + Mfl £ f ^ (A1Q) 
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with 



Pa 



r(Yo) 



dfl a 



(All) 



In order to proceed beyond the "no Coulomb" approxi- 
mation one should perform the expansion in 1A2|I at least 
including the quadratic terms. With these terms one ar- 
rives at the linearized Poisson equation for the Coulomb 
potential that incorporates the Debyc screening, cf. H0. 



is non-zero. Thus the averaging region for proton and 
kaon is inside the lump and that for electron is the cell. 

To obtain correct values of the Debye screening 
lengths, the implicit dependence of (p a ) on p a through 
the change of mean fields and the Coulomb potential 
should be taken into account. We will see that it is es- 
pecially in the case of the kaon Debye term. Indeed, the 
linearization procedure should be performed in all fields 
and the linearized Poisson equation is as follows: 



APPENDIX B: DEBYE SCREENING LENGTH 

To roughly estimate the screening effect we can intro- 
duce the Debye screening lengths as follows. A linearized 
Poisson equation for the Coulomb potential (SVcoui = 
Vcoui — fi>e ~ Vq, where Vq is the reference constant value, 
takes the form |70| 

A^coui + (A^ ot) ) ~ 2 SVcoui = 0, (Bl) 

(Ag* J )" a = E ^ 

a—e,p,K 

A («)J te2 «" /2 (B3) 

V dfJ-a J 

The partial contributions A^' , the Debye screening 
lengths, are dependent on the region. To find typical val- 
ues characterized screening inside the structure we aver- 
age p a over the region where the relevant charge density 



A^Vcoui + avvSVcoui + a Vrj 5a (B4) 
+ay u fo + avpSRo + a V g59 
+avn5p n + avpSpp = 0, 

where ayi are density dependent coefficients. The fields 
5a, Sloq, 5Rq, 59 are determined by their own equations 
of motion having similar forms. Also variation of equa- 
tions for p n and p p yields equations for Sp n and 5p p , p e 
is unambiguously expressed via the neutron and proton 
chemical potentials, therefore we did not present the cor- 
responding term in (|B4fl . Neglecting derivative terms for 
all the fields except for the Coulomb field in the system of 
equations similar to (|B4fl one can express 5a, Sioq, 5Ro, 
Sp n and Sp p as functions of the only one SVcoui variable 
and one finally arrives at Eq. 1|B1|) . However such a pro- 
cedure is legitimate only for short-scale fields, 8a, Sujq, 
SRq, Sp n and Sp p but not for the kaon field 89 since the 
later field proves to be a long-scale field. The value AS9 
entering the kaon equation of motion cannot be omitted 
thereby. 
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